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1. Introduction 

The existence of the standard model Higgs boson is a necessary ingredient for the consistency 
of the standard model. Although the mass of the Higgs boson cannot be predicted, bounds on the 
Higgs boson mass can be given. A lower bound can be derived from the requirement of a stable 
vacuum while an upper bound can be computed from the triviality of the theory which leads to the 
occurrence of the Landau pole. These mass bounds are usually computed in perturbation theory and 
indeed, perturbative results for the Higgs boson mass bounds were obtained in the past. However, 
the perturbative analysis may be questioned, since for the upper bound, the quartic coupling can 
become strong. The lower bound, on the other hand, may be an artefact of perturbation theory 
when looking at Higgs field values far away from the minimum of the potential. These concerns 
motivated a non-perturbative and ab initio lattice field theory investigation of the mass bounds of 
the Higgs boson in a chirally invariant lattice Higgs- Yukawa model [|l|, ^. 

In the lattice approach it is straightforwai^d to detennine those mass bounds also for the case 
that there is a fourth, heavy family of quarks. Such an extension is not excluded a priori and further, 
a fourth family of fermions offers the possibility of generating sufficient CP violation to fulfil the 
Sakhai^ov condition to explain electroweak baiyogenesis of the early universe [Q]. 

However, CKM4 fits and direct searches already make a naive fourth generation scenario quite 
improbable [Q], especially if the newly discovered particle with a mass around 126 GeV ^ turns 
out to be the standard model Higgs boson. In this work we provide additional general constraints 
on the fourth family quark masses by their effect on the Higgs boson mass bounds. 

2. Model and Implementation 

We study the Higgs- Yukawa sector of the standai^d model. The field content we consider is a 
fermion doublets *F and a complex scalar doublet <p. The continuum action in this model is given 
by: 

+ J d'^x{t^t + b^b + ytoYL9bi,+yto¥L(pt^ + h.c.} , with ^ = /T2<p*. (2.1) 

The bare quartic self coupling of the scalar field is given by Xq, niQ denotes the bare mass of the 
scalar field and ytg/ho the Yukawa couplings of the fermion fields. We want to stress, that we do 
not include any gauge fields in this model neither the gluonic degrees of freedom nor electroweak 
ones. This is done for computational simplicity and further we expect the effect of the gauge fields 
to be small for the problems addressed in our work. 

We will mostly study the system in the phase with spontaneous symmetry breaking where the 
scalar field develops a non vanishing vacuum expectation value (vev) and three Goldstone modes 
emerge. Further we only consider the heaviest fermion doublet, since the dynamics of the scalar 
field is dominated by the largest Yukawa coupling. 
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For the discretized lattice version of the bosonic action we rewrite the scalar doublet as a real 
four vector <I>. With the repaiametrisations 

/— + , X 2 2 I-2X-SK 

a(p = V2K\„_.^^\, ^=4^' «"Jo = (2-2) 

the bosonic action can then be written in the compact form: 

= -kI^^I [<I>,+^ +0,-^] +£<I>J<I>, + A£ [OtO,- 1]'. (2.3) 

X,II X X 

Here K denotes the hopping parameter. The subscripts x and x^pL on the scalar fields denote the 
scalar field at the space-time point x and x±afi respectively. 

For the discretization of the fermions we use the overlap operator D"^' [0, |8|, ^ with a Wilson 
Dirac kernel : 

Z)- = p{l + -A=|, A = D^-p, D^ = l^r,y',-'^y',yi, (2.4) 



with V;^'^"^ being the forward, backward and symmetrized lattice nearest neighbor difference op- 
erators in direction jj.. The so-called Wilson parameter r is as usual chosen to be r = 1. The 
dimensionless parameter p is free to be chosen in the range < p < 2r. However, the locality 
properties of the free overlap operator in case of vanishing gauge couplings are optimal for the case 
of p = 1 which will therefore be the choice in our work. The fermionic part of the action is 
then: 

Sf[y, Y,^]=J^ ^x [D''" + PJ^diag{y„yb)P+ + P_diag{y,,yb)(t>xP^] Vx, St/h = VlKyt/h,- 

(2.5) 

This action now obeys an exact global SU(2)l x U(l)y lattice chiral symmetry with the transfor- 
mations: 

Y^UyP+Y + Uy^lP-V, W^W+^l^y + W-Uy, (p^UYip^l, ^+ ^ HL^ff/^, (2.6) 
for any D.i G SU(2)/, and Uy € U(l)y. The modified chiral projectors are given by: 

P,_ = i±^, f = ^(l-^D-) (2.7) 



2 ' ■ • V P / 

Even though in principle different masses for the fermions in the doublet are possible we 
restrict ourselves to a mass-degenerate doublet in this work. For the implementation we use a 
polynomial Hybrid Monte Carlo algorithm [10] with various improvements implemented. For 
further details of the implementation see [11]. 

To set the scale, i.e. to determine the lattice spacing a, we use the phenomenologically known 
value of the vev of 246 GeV. Further we define the cutoff A as the inverse of the lattice spacing. 
Since in finite volume a naive average scalar field would vanish without an external field, we define 
the magnetization m as the average absolute value of the scalar field: 



m 



\^\), ^ = V 'Y.^,, (2.8) 
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with V being the volume of the space time lattice. This approach has the same thermodynamical 



limit as the application of an external source []12|]. The magnetization and the renormalized vev are 
related as follows: 



m\/2K 



a 



246 GeV, A : 



246 GeV 



(2.9) 



The Goldstone and Higgs bosom renormalization constants Z^/^ and their masses are computed 
from the real part of the Goldstone or Higgs boson propagators Gq/h respectively: 



0. 



(2.10) 



The propagators are computed for discrete lattice momenta and fitted according a one loop moti- 
vated formula derived in lattice perturbation theory. It may be noted, that this approach especially 
for the Higgs boson is only valid if the decay width of the particle is small compared to its mass. 
However, it was shown in a rigorous resonance analysis [|3|, that at least in the case of a physical 
top quark mass, the pole of the propagator fully agrees with the resonance mass determined in [p^]. 
Finally, the masses of the fermions nif are computed by means of the time slice correlator 



Cf{At) 



L,-l 



£ l^{23iTr{PMt + ^t,x)-Y{t,y)P,}), 



(2.11) 



which shows a behavior proportional to cosh{amf{At — L,/2)) for large time separations At. 

In general we observe rather severe finite volume effects due to the almost massless Goldstone 



modes, which cause finite volume effects proportional to Lj^ [14, 15| rather than an exponential 
falloff with Euclidean time for theories with a mass gap. Some examples for finite volume effects 
can be found in fig. |l} There, one can see, that for a trustworthy determination of the Higgs boson 
mass and the vev lattice sizes of up to at least Ls = 32 are necessary. 

In addition to the non-perturbative determination we also perform a perturbative analysis by 



means of the constrained effective potential as described in detail in [11]. Those calculations are 
perfomied within the same lattice regularization as the numerical simulations. We employ discrete 
lattice sums for the loop corrections which are then evaluated numerically. The effective potential 
V{^) is determined to one loop in the large Nf limit with the vev and the Higgs boson mass given 
by: 



0; 



m 



H- 



(2.12) 



In order to have a stable vacuum in the scaling regime (|0| < 0.5), we demand the potential to be 
concave up everywhere in this region. In this perturbative framework we can determine the lower 
Higgs boson mass bound by finding the lowest possible Higgs boson mass while keeping physical 
quantities fixed and still fulfil the required stability conditions. 



3. Results 

Let us shortly summarize the strategy to determine mass bounds for the Higgs boson in our 
analysis. In total we have three free parameters, namely the bare quartic self coupling X, the 
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(b) top quark mass 
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(c) Higgs boson mass 



Figure 1 : We show finite volume effects for the magnetization, the quark mass and the Higgs boson mass at a 
cutoff of about 1 .5 TeV. The data shown correspond to simulations performed at infinite bare quartic coupling 
and fermion masses in the range m,/ « 200. . .700 GeV. The lattice sizes shown are = 12, 16,20,24,32. 
The red lines indicate linear fits involving only data from computations with > 20 the grey lines show 
quadratic fits over the whole range of lattice extents. 



hopping parameter K and the Yukawa coupling y. We fix the Yukawa coupling such to obtain the 
desired mass for the fermions and tune K to retrieve the cutoff at which we want to investigate the 
mass bounds. Then it can be shown [[T], that with the cutoff and the fermion mass fixed, the 
smallest accessible Higgs boson mass is obtained for choosing A = while the Higgs boson mass 
is the largest for X =o°.ln earlier works we investigated the cutoff dependence of the Higgs boson 
mass bounds for the case of a fermion doublet at the physical top quark mass rrit 175 GeV ^ 



and for a very heavy doublet with a mass around nifi 676 GeV [|16p. The results can be found in 
fig. While the upper Higgs boson mass bound is only increased slightly, the lower mass bound 
for the Higgs boson increased by a factor around 5 when the fermion mass is increased. To obtain 
a better understanding of the increase in the mass bounds, we investigated those bounds for several 



fermion masses at a cutoff around A 1.5 TeV. Those results are shown in fig. ^ One clearly 
observes the smooth increase of the lower mass bound of the Higgs boson with increasing quark 
mass. Our results suggest that a Higgs boson mass of 126 GeV would allow heavy quarks only up 
to approximately 300 GeV. 

In fig. ^ we also show the lower Higgs boson mass bound obtained from the perturbative 
effective potential calculations. It is remarkable, how well the perturbative result agrees quali- 
tatively with the non-perturbative findings up to quark masses of about 700 GeV (corresponding 
to a Yukawa coupling around 3). This lets us assume, that the model may be perturbative in a 
wide region of the parameter space. This finding allows us to test perturbatively whether adding 
a dimension-6 operator, i.e. a Xe,(p^ term, to the bosonic action may have an impact on the lower 
bound. To this end, we determined the lower bound at a cutoff A = 2 TeV for various quark masses 
with and without a dimension-6 operator at two different couplings Ae. Those results are shown in 
fig. ^ We clearly see that the additional operator has no visible impact on the lower Higgs boson 
mass bound. 
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(a) physical top quark mass 
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(b) heavy quark double 



Figure 2: We show the cutoff dependence of the Higgs boson mass bounds for two masses of a mass 
degenerate quark doublet. On the left graph we show the data for the physical top quark mass while the right 
graph corresponds to a very heavy mass around 676 GeV. 
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(a) Results at A ~ 1.5 TeV (b) Higher dimensional operators 

Figure 3: We show Higgs boson mass bounds at fixed cutoff. On the left graph non-perturbative data are 
shown for the upper and lower mass bound. (For the red triangles the infinite volume extrapolation has 
not been carried out yet). The data points for the lower bound are compared to findings in a perturbative 
effective potential calculation. The right-hand plot compares perturbative results for the lower bound with 
and without a dimension-6 operator included at a cutoff A = 2 TeV. 



4. Conclusion and Outlook 

We performed a non-perturbative determination of Higgs boson mass bounds in a Higgs- 
Yukawa model with a physical top quark mass and a heavy quark doublet. We found that a 126 GeV 
Standard Model Higgs boson excludes a naive fourth generation of quai^ks if their mass exceeds 
300 GeV. Further we gave perturbative arguments, that this bound is stable against the inclusion of 
higher dimensional operators. However, we plan to addressed this finding also non-perturbatively 
in the future. Further, our setup allows to test whether a non -degenerate doublet could alter the 



mass bounds, since a mass splitting of up to mi//mt' « 0.9 is possible |17]. Another interesting 
direction would be to include an additional scalar field, as suggested in [18], which could have an 
impact on the Higgs boson mass bound. 
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